APPENDIX C 

/* Normalized Adaptive Noise Canceler */ 



#include <stdlab.h> 
#include <math.h> 



#define TRUE 
#define FALSE 

#define FLOAT 3 2 
#define INT32 
#define VOID 

#define FABSF 
#define SQRTF 

/define MAX(a,b) 
/define MIN(a,b) 



1 
0 

float 
long int 
void 

f abs 
sqrt 

(a) > (b) ? (a) : (b) 
(a) < (b) ? (a) : (b) 



/define MIN_VAL 
/define MAX_DEL 
/define MIN_DEL 
/define MAXRHO 
/define MIN_RHO 
/define MIN BSERR 



0.01 

0.999999 
-0.999999 

2.0 
-2.0 

1E-15 



typedef struct { 



FLOAT3 2 


berr ; 




FLO AT 3 2 


berr_l 




FLO AT 3 2 


delta ; 




FLOAT 3 2 


err; 




FLOAT3 2 


f err ; 




FLOAT3 2 


gamma; 




FLO AT 3 2 


gamma_ 


i; 


FLO AT 3 2 


rho; 




FLO AT 3 2 


delta_ 


i; 


FLOAT 3 2 


Bserr ; 




FLO AT 3 2 


Bserr__ 


i; 


LANC CELLS; 







typedef struct { 

INT32 cc; 

FLOAT3 2 lambda; 

FLO AT 3 2 min_error; 

LANC_CELLS *cells; 
} LANC Context; 



/* number of cells 

/* put in value for lambda 

/* parameter 

/* point to array of ANC__CELLS 



extern LANC_Context * 

LANC__Init( 

INT3 2 num^cells, 
FLO AT 3 2 lambda, 
FLOAT 3 2 min error) ; 



/* number of cells 
/* lambda param 
/* min error 



extern VOID 
LANC_Done ( 

LANC Context 



*c) ; 



extern VOID 
LANC_Reset ( 

LANC Context 



*anc) ; 



extern FLO AT 3 2 
LANC_Calc( 

LANC_Context *anc, /* input, context handle */ 

FLO AT 3 2 rips, /* input, noise plus signal */ 

FLOAT32 noise) ; /* input, noise reference */ 

/* The following macros provide efficient access to the lattice */ 



#def ine 


ANC_CELL 


_SIZE 


11 








#def ine 


xBERR 


0 










/define 


xBERR_l 


1 










#def ine 


xDELTA 


2 










#def ine 


xDELTA 1 


3 










/define 


xGAMMA 


4 










/define 


xGAMMA_l 


5 










/define 


XBSERR 


6 










/define 


XBSERR_1 


7 










/define 


XERR 


8 










/define 


XFERR 


9 










/define 


xRho 


10 










/define 


berr 
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/define 
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/define 


P_berr 
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#def ine 


berr 1 




( it 
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/define 


Bserr 
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/define 


Bserr_l 
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(P 
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/define 


P_Bserr_ 
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/define 


ir delta 
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/define 


err 
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/define 


P_f err 




(* 


(P 


+ 


XFERR - ANC_CELL_SIZE) ) 


/define 


f err 




(* 


(P 


+ 


XFERR) ) 


/define 


gamma 




(* 


(P 


+ 


xGAMMA) ) 


/define 


P_gamma 




(* 


(P 


+ 


XGAMMA - ANC CELL SIZE)) 


/define 


N_gamma 




(* 


(P 


+ 


xGAMMA + ANC CELL SIZE)) 


/define 


P_gamma_ 


1 


(* 


(P 


+ 


XGAMMA 1 - ANC CELL SIZE)) 


/define 


gamma_l 




(* 


(P 


+ 


XGAMMA1) ) 


/define 


rho 




(* 


(P 


+ 


XRho) ) 



/* 

Name: LANC_Init 

Abstract: Create an ANC context 



extern LANC Context * 
LANC lnit( ~ 

INT32 num_cells, /* number of cells */ 

FLOAT 3 2 lambda, /* lambda par am */ 

FLOAT3 2 min_error) { /* min error */ 

LANC_Context *anc; /* context */ 

anc = (LANC_Context *)malloc(sizeof (LANC Context) ) ; 
assert(anc != NULL) ; 
anc->cc = n\im_cells; 
anc->lambda ~ lambda; 
anc->min error = min error; 

anc->celTs = (LANC_CELLS * ) malloc (sizeof (LANC_CELLS) * (num_cells + 2)); 
assert (anc->cells != NULL); 

return (anc) ; 



/* 



Name: LANC_Reset 

Abstract: Reset an ANC context 

*/ 



extern VOID 
LANC_Reset ( 

LANC_Context *anc) { 

FLO AT 3 2 *p; 
INT3 2 m; 



p = ( FLO AT 3 2 *)anc->cells; 

for (m = 0; m <= anc->cc; m++) { 



rho 






o . 0 ; 




err 






0.0; 




f err 






0.0; 




berr 






0.0; 




berr_l 






0.0; 




delta 






0.0; 




delta_ 


1 




0.0; 




Bserr 






anc->min_ 


error ; 


Bserr_ 


1 




anc->min 


terror ; 


gamma 






MIN VAL;~ 




gamma_ 


1 




MIN VAL; 




p 




+= 


ANC CELL 


SIZE; 



p = (FLOAT32 *)anc->cells; /* Cell # 0 special case */ 

gamma = 1.0; 
gamma_l = 1.0; 

} 



/* 

Name : LANC_Done 

Abstract: Delete an ANC context 



*/ 



extern VOID 
LANC_Done ( 

LANC Context 



*anc) { 



free(anc->cells) ; 
free(anc) ; 



} 



/* 



Name: LANC Calc 



Abstract : Calculate 



FLOAT 3 2 



LANC_Calc( 

LANC_Context *anc, /* input, context handle 

FLOAT32 rips, /* input, noise plus signal 

FLOAT32 noise) /* input, noise reference 



*/ 
*/ 
*/ 



{ 

INT32 m; 

FLO AT 3 2 *p; 

FLOAT32 B,F,B2,F2; 

FLO AT 3 2 qd2,qd3; 

I NT 3 2 output_cel 1 ; 



/* Update time delay elements in cell structure 



*/ 



p = (FLO AT 3 2 *) anc->cells; 

for (m = 0; m <= anc->cc; m++) { 

gamma_l = gamma; 

berr_l = berr; 

Bserr__l = Bserr; 

delta_l = delta; 

p += ANC_CELL_SIZE; 

} 

/* Handle Cell # 0 */ 

p = (FLO AT 3 2 *)anc->cells; 

Bserr = anc->lambda * Bserr_l + noise * noise; 
Bserr = MAX (Bserr, MIN_BSERR) ; 

ferr = noise / SQRTF (Bserr) ; 
ferr = MAX (ferr, MIN_DEL) ; 
ferr = MIN(ferr, MAX_DEL) ; 

berr = ferr; 

rho = anc->lambda * SQRTF (Bserr_l / Bserr) * rho + berr * nps; 
N_err = nps - rho * berr; 

output_cell = anc->cc - 1; /* Assume last cell for starter */ 

for (m = 1; m < anc->cc; m++) { 
p += ANC_CELL_SIZE; 

B = SQRTF (1.0 - P_berr_l * P_berr_l) ; B2 = 1.0/B; 



F = SQRTF (1.0 - P^ferr * P_ferr ); 



F2 = 1.0/F; 



P_delta = P_delta_l * F * B + P_berr 1 * P_ferr; 
P_delta = MAX(P_delta, MIN_DEL) ; 
P_delta = MIN(P_delta, MAX_DEL) ; 



qd3 
qd2 

f err 
f err 
f err 

berr 
berr 
berr 



= 1-0 - P_delta * P_delta; 
= 1.0 / SQRTF (qd3 ) ; 

= (P_ferr - P_delta * P_berr_l) * qd2 * B2 ; 
= MAX (f err, MIN_DEL) ; 
= MIN(ferr, MAX_DEL) ; 



= (P_berr_l - P_delta * PJerr 
= MAX(berr, MIN_DEL) ; 
= MIN(berr, MAX_DEL) ; 



) * qd2 * F2; 



} 



gamma = P_gamma * (1.0 - P_berr * P_berr) ; 
gamma = MAX (gamma , MIN_VAL) ; 
gamma = MIN (gamma, MAX_DEL) ; 

Bserr = P_Bserr_l * qd3; 

Bserr = MAX (Bserr, MINBSERR) ; 

rho *= anc->lambda * SQRTF ( (Bserrl / Bserr) * (gamma / gamma_l) ) ; 

rho += berr * err; 

rho = MAX(rho, MIN_RHO) ; 

rho = MIN(rho, MAX_RH0) ; 

N_err = err - rho * berr; 



p = (FLOAT32 *) &(anc->cells[output_cell /* *ANC_CELL SIZE */]); 
return (N_err) ; ~~ 



